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We describe a simple mechanical system, a ball rolling along a specially-designed landscape, that mimics 
the dynamics of a well known phenomenon, the two-bounce resonance of solitary wave collisions, that has 
been seen in countless numerical simulations but never in the laboratory. We provide a brief history of the 
solitary wave problem, stressing the fundamental role collective-coordinate models played in understanding 
this phenomenon. We derive the equations governing the motion of a point particle confined to such a surface 
and then design a surface on which to roll the ball, such that its motion will evolve under the same equations 
that approximately govern solitary wave collisions. We report on physical experiments, carried out in an 
undergraduate applied mathematics course, that seem to verify one aspect of chaotic scattering, the so-called 
two-bounce resonance. 

PACS numbers: 05.45Yv, 05.45Gg 

Keywords: solitons, chaos, solitary wave collisions, mechanics 


Solitary waves are solutions to partial differen¬ 
tial equations that maintain their spatial profile 
while moving at constant speed. A wide variety of 
systems display a behavior called chaotic scatter¬ 
ing when two such waves collide. The two waves 
may bounce off each other one or more times be¬ 
fore escaping to infinity, or they may capture each 
other and never escape. The number of colli¬ 
sions and the final speed of those that separate 
depends in a complex way on their initial speeds. 
To our knowledge, this process has never been ob¬ 
served in a laboratory experiment involving soli¬ 
tary waves. The problem is described well by a 
small finite-dimensional system of ordinary dif¬ 
ferential equations. We describe an experiment, 
a ball rolling on a specially designed surface, that 
obeys the same system of ODE. We report on 
laboratory experiments demonstrating that the 
experimental system has similar dynamics to the 
solitary wave collisions. 


I. INTRODUCTION 

Solitary waves , localized structures that translate at 
constant velocity while maintaining their spatial pro¬ 
file, are among the most important concept in nonlin¬ 
ear physics. Solitons are a class of solitary waves which 
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possess an additional underlying mathematical structure, 
one which renders the equations solvable, and which, in 
particular, causes the behavior of colliding solitons to be 
strikingly simple: solitons survive a collision with their 
shape and velocity intact, but with their positions shifted 
by a finite, computable, amount. 

As long as scientists have known about solitary waves 
and had computers capable of simulating them, we have 
been colliding solitary waves together numerically and 
observing what comes out. An important example is the 
collision of so-called kink and antikink solutions to the 
equation, 

fa - Vxx + - T = 0- (i) 

The first studies in the 1970’s gave a fleeting hint of 
rich structure and the first step toward understanding 
it. This was explored more deeply and given a name, 
the two-bounce resonance , in the 1980’s. In the 1990’s, 
more detailed numerical experiments revealed chaotic 
scattering. Finally, in the early 2000’s, the mechanism 
behind this chaotic scattering was explained fully using 
techniques from dynamical systems. Figure [l] shows the 
speed at which a kink and antikink separate as a func¬ 
tion their speed of approach, demonstrating a rich struc¬ 
ture. Numerical simulations of partial differential equa¬ 
tions (PDE) arising in diverse areas of physics have re¬ 
vealed the same phenomenon. In all that time, however, 
this behavior has never been reported in a physical exper¬ 
iment in a real nonlinear wave system. 

This paper describes a simple experiment—a ball set 
rolling on a manufactured landscape—which mimics the 
dynamics of solitary wave collisions and reproduces some 
features of Figure [l] The main tool for understanding 
solitary wave interaction has been the derivation and 
analysis of collective coordinate (CC) models. The land- 
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FIG. 1. The escape speed as a function of the input speed 
for a kink-antikink collision in the ip 4 system for PDE |T), 
phenomenological ODE system and the discrete map ( fl2| ) . 
Color indicates number of collisions before separation: one 
(black), two (blue), three (green), and four (red). 

scape is constructed such that the ball’s equations of mo¬ 
tion strongly resemble the ordinary differential equation 
system (ODE) governing the evolution of solution pa¬ 
rameters described by a CC reduction describing solitary 
wave collisions in the equation. 

Chaotic scattering is an appealing phenomenon to 
study because images such as Figure [l] cry out for ex¬ 
planation. The mathematical structure that has been 
developed is equally appealing, as is well-described by 
OttP 2 A scattering process is a physical phenomenon in 
which an object’s trajectory begins and ends in free mo¬ 
tion (i.e. with zero acceleration) but spends a finite time 
in a region where it is subject to forces. Such a process 
is called chaotic scattering if the final state depends in a 
complex, i.e. fractal or multi-fractal, manner on its initial 
state. 

The paper is organized as follows. Section [TT] discusses 
where the problem comes from and what approaches have 
been most useful in understanding it. It begins with a 
brief recap of the history of solitons, followed by an ex¬ 
planation of the difference between solitons and garden- 
variety solitary waves. It then describes the long history 
of chaotic scattering of solitary waves, the progress that 
has been made, and that methods that have been used 
to explain the observations. In Section m we describe 
the laboratory experiment and the mathematical model 
that describes it, while in Section |TV| we present the ex¬ 
perimental results, which we interpret as the two-bounce 
resonance, the aspect of chaotic scattering that is most 
robust in the presence of dissipation. We end in Section|V| 
by reviewing the physical systems that have motivated 


many of earlier studies of the two-bounce/chaotic scatter¬ 
ing phenomenon and expressing our hope that someone 
may perform this experiment to display this phenomenon 
in an extended physical system. 


II. HISTORICAL MOTIVATION 
Solitary waves and solitons 

Solitary waves are solutions to a nonlinear evolution¬ 
ary partial differential equation (PDE) that translate at 
a constant speed while maintaining its spatial profile. 
Solitons are solitary waves which exhibit unexpectedly 
simple dynamics upon collision. This is due to a hidden 
mathematical structure that was discovered beginning in 
the 1960’s, first for the Korteweg-de Vries equation®!, 

u t + 6uu x + 

^XXX — 0 } 

and later for many other equations, such as the sine- 
Gordon equation, 

u u ~ u xx + sin u = 0. (2) 

Such soliton equations can be solved exactly using a 
method called the inverse scattering transform (IST)P 
Good brief histories of solitary waves and solitons can be 
found elsewhere, e.g. in encyclopedia articles by ScottP, 
and by Zabusky and Porter! 7 

The 1ST can be used to explicitly show that after col¬ 
liding, solitons continue propagating with the same shape 
and speed, but with a phase shift and a time delay. Most 
nonlinear wave equations do not possess the 1ST, how¬ 
ever, and solitary waves collisions in such systems display 
more complicated dynamicsn. Beginning in the 1970’s, 
numerical simulations hinted this picture, but it took the 
development of more powerful computer hardware and 
numerical methods for the true complexity of this dy¬ 
namics to fully emerge. The simultaneous development 
of collective coordinate methods would prove useful for 
untangling this complexity. 


Collective coordinates methods 

The term “collective coordinates” refers to a variety 
of methods that are used to obtain simple ODE mod¬ 
els that approximate the behavior of a PDE, at least as 
long as the solution stays in some small region of solu¬ 
tion space. A very useful approach is the “variational 
method”, which applies to partial differential equations 
that arise as the Euler-Lagrange equations for a system 
with Lagrangian density 

I = C((p,x,t)dx dt (3) 






























































3 


For example, solutions of the Lp A equation 0 minimize 
the action due to 

r* 1 2 1 2 , 1 2 1 4 

£= 2 ^“ 2 ^ + 2 ^ - 4 ^ 

over all C 1 function satisfying appropriate boundary con¬ 
ditions at infinity. 

The variational method, due originally to Bondesonf^ 
works by constructing a solution ansatz whose spatial 
profile depends on a small number of time-dependent 
parameters and minimizing the action § with respect 
to these parameters. The resulting equation is a finite¬ 
dimensional Lagrangian system of ODE governing the 
parameters’ evolution. Many examples and some gener¬ 
alizations are given in the review paper of Malomed. 9 


Numerical and analytical studies of solitary wave collisions 


Solitary waves exist in the (/? 4 equation ([!]), which is 
non-integrable, despite its clear similarity to the sine- 
Gordon equation These “kink” waves are given by 


ifK (x, t', v) = tanh —=, where £ = 
v 2 


x — vt 
y/l — V 2 


for any — 1 < v < 1. (4) 

The “antikink” solution is just cp^(x,t,v) = —<^k {x,t\v). 

Many papers have been written about kink-antikink 
(KK) collisions in this system. In 1975, Kudryavtse\ff^ 
found numerically that when v = 0.1, the KK pair co¬ 
alesce into a single bound state at the origin which os¬ 
cillates irregularly and is slowly damped to zero. He ex¬ 
plained this with a formal argument based on potential 
energy which can be considered a first step toward an 
explanatory collective coordinates model. 

In 1979, Sugiyamsf 11 performed additional experiments 
for a few values of v between 0.1 and 0.6. He found that 
those with speeds below a critical value v c « 0.25 the 
KK pair merge into a localized bound state, while for 
v > u c , the pair collides inelastically. Careful examina¬ 
tion of his simulations showed that, upon collision, an os¬ 
cillatory mode (sometimes called a “shape mode”, since 
its effect is to alter the shape of the kink) is excited, 
which removes energy from the translation component. 
If the initial kinetic energy is less than the amount lost 
to the oscillatory mode, the KK pair can not escape to 
infinity. If the v > u c , they will escape to infinity but 
with reduced speed and with oscillations superimposed. 
He also derived the first CC model for this phenomenon, 
and used it in a formal calculation to determine v c in 
agreement with his numerical observations. 

The shape mode is an eigenfunction for the lineariza¬ 
tion of equation ([l]) about the kink solution Q, and is 
given by 


xi(0 


1 





FIG. 2. (a) The potential and coupling function of equa¬ 

tion (b) The phase plane of the X — X subsystem, sepa- 
ratrix given by thick line. 


and oscillates with frequency uj\ = y Sugiyama based 
his CC ODE for this system on the ansatz 


tp(x,t) =lf K 0 - X(t)) + If K ( x + x(t)) - 1 

+ A (t)x i (x - x(t)) - A(t)x 1 (x + X(t)). 


( 5 ) 


It leads to a two degree-of-freedom Hamiltonian system 
for the evolution of X(t) and A(t). The ODE system that 
arises is somewhat complicatedf- 11 but its fundamental 
features are captured by the simpler model system 

X(t) + U\X ) + cAF'(X) = 0; (6a) 

A(t)-hcu 2 A-hcF(X) = 0, (6b) 

where U(X) = e~ 2X — e~ x and F(X) = e~ x , plotted 
in Figure (2|a). This conserves an energy 

E=^X 2 + U(X) + I (i 2 + w 2 A 2 ) + cAF(X). 


The term U(X) is a potential energy describing the in¬ 
teraction of the two kinks, and was essentially derived 
by KudryavtsevP For large X, both F'(X) and U'{X) 
vanish, so that the kink and antikink centers XX (t) move 
at constant speed. The most important feature, for our 
purposes, is the phase diagram for X when c = 0, shown 
in Figure |2|b). The trajectories are level sets of the en¬ 
ergy E =!X 2 + U(X), with E < 0 on bounded orbits, 
a separatrix orbit with E — 0 and unbounded orbits for 
E > 0. We denote the separatrix orbit by X$(t), and 
note it is an even function of time. 

Figure [l] shows that this model qualitatively captures 
much of the dynamics of solitary wave collisions in ([!). 
Observe, however, a few fundamental differences. While 
both systems conserve energy, in extended system ([!), 
radiation (phonons) can carry energy away from the im¬ 
mediate vicinity of the kinks, effectively adding dissipa¬ 
tion. This can be seen in the figure, where only solutions 
that escape after four or fewer collisions are plotted. In 
ODE Q, all solutions escape except for a set with mea¬ 
sure zero, while for the PDE, a significant fraction of so¬ 
lutions are trapped, due to energy loss to radiation. The 
radiation also steals some energy from solutions that do 
escape, so that v out < V[ n even at the maxima of the 
resonance windows. 
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Also in 1979, Ablowitz, Kruskal, and Ladik reported 
on a series of numerical experiments for a few differ¬ 
ent nonlinear Klein-Gordon equations of the form ip tt — 
<p xx + fi}P) = 0, including the cp 4 equation and the sine- 
Gordon equation. 12 . Their numerical results were similar 
to Sugiyama’s but the paper ends with an observation 
which was to prove important: On certain velocity inter¬ 
vals below f c , the kink and antikink eventually separate, 
instead of forming a bound state as seen by Kudryavt¬ 
sev (from their References, it seems they had not seen 
Sugiyama’s work at this point). When the initial veloc¬ 
ity is 0.3—greater than v c —the final velocity is 0.135. 
However, when the initial velocity is 0.2—less than v c , 
the final velocity is 0.155. They conclude “The reason 
for the apparent ‘resonance’ between these interacting 
aperiodic waves and the radiation is not yet fully under¬ 
stood.” 

In 1983, Campbell et al. 13 performed a more system¬ 
atic numerical sweep of initial velocities and found sig¬ 
nificantly more detailed structure. Their calculation re¬ 
vealed the black curve and the leftmost nine blue curves 
in Figure^ The black curve shows the final velocities of 
all the solutions that escape after exactly one collision, 
i.e. those with initial speed above v c , which they esti¬ 
mate numerically to be 0.26. The blue curve shows the 
final velocities of those with exactly two collisions. They 
called this phenomenon the two-bounce resonance. 

The first question is what is the difference between 
solutions in one two-bounce window and those in the 
next. A good way to understand the solutions is to 
fit the numerical solution to the ansatz <R and plot 
the approximate values of X(t) and A(t). This is done 
in Figure [3] for five increasing initial velocities. Subfig¬ 
ure (a) shows a collision leading to capture; subfigures (b) 
and (d) show solutions from the first two two-bounce win¬ 
dows; (c) shows a solution from a three-bounce window; 
(e) shows a one bounce solution with V[ n > v c . The 
number of oscillations of A(t) is the same for all initial 
velocities in a given window and increases by one from 
one window to the next. Let T s h ap e be the period of the 
shape oscillation and T n be the time interval between the 
two collisions for the v- m in the nth window at which v out 
is maximized, with T 2 and T 3 indicated in subfigures (b) 
and (d). We remark that the nonlinear projection (J5| is 
singular when X(t) = 0, because the coefficient of A(t) 
then vanishes pl in which case the fitting algorithm fails. 

Treating the numerical results as laboratory data, 
Campbell et al. reasoned there should be a relation of 
the form 


T n ~ T s h a p e n T 5) 


( 7 ) 







FIG. 3. Position X(t) (solid, left y- axis) and shape-mode 
amplitude A(t) (dashed, right y- axis) for selected initial ve¬ 
locities. (a) v = 0.184, Capture, (b) v — 0.1986, First two- 
bounce window, (c) vq = 0.2236, a three-bounce window, 
(d) v = 0.2268, second two-bounce window, (e) v = 0.27 
escape without capture. Vertical asymptote in (a) due to sin¬ 
gularity of finite dimensional reduction, as per Caput 


detailed experiments, the existence of higher bounce win¬ 
dows arranged in a fractal structure, i.e. of chaotic scat¬ 
tering, in the solitary wave collisions, including now, the 
three and four bounce windows of Figure [l] 

Over the next two decades, these numerical experi¬ 
ments were reproduced for many other nonlinear wave 
systems, and Campbell et al.’s reasoning was applied to 
these systems as well. Campbell, as well as Peyrard, 
looked at a variety of (p 4 -like equations. They demon¬ 
strated that, at least for the systems they consid¬ 
ered, the resonance phenomenon was only present when 
the linearization around the kink possesses an internal 
mode. 16 18 Of particular interest is a study by Fei et all 19 
which considers the sine-Gordon equation perturbed by 
a small localized defect, 


which they confirm by least-squares fitting the data to a 
line. They conjectured that at the first collision, the kink 
loses energy to the shape-mode oscillation, and that, if 
the timing is right on the second collision, the shape mode 
returns enough energy to the propagating mode to allow 
escape. What remained was to explain the mechanism. 

A few years later, Anninos et alP^ showed, via more 


Utt ~ u xx + sinu m eS(x) sinu. (8) 

They find that kink solutions impinging on the defect 
show a behavior very similar to that seen in Figure [l] 
They derive a CC system similar to system © where 
A(t) now measures the amplitude of a mode localized in 
neighborhood of the defect. They use conservation of 
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energy to derive an implicit formula for v c (e) and ap¬ 
ply Campbell et al.’s reasoning to fit the time between 
collisions to the period of the secondary oscillator. 


Reduction of collective coordinates to an iterated map 

Beginning in 2004, Goodman and Haberman published 
a series of papers explaining the chaotic scattering in 
detail by reducing the CC ODE system to a discrete¬ 
time Iterated map, called a separatrix map or scattering 
mapPS 21 This was first done^ for the system studied 
by Feip^ because the the CC ODE system for equa¬ 
tion © has a small parameter e that can be used in 
a perturbation analysis. Explicit formulas were found 
to approximate the critical velocities and the two- and 
three-bounce resonance, eliminating the need for data 
fitting as in equation (|7|). Over subsequent papers pMUl 
the derivation of the map was streamlined and applied to 
other solitary waves systems and the model system li¬ 
lt was subsequently put in a more explicit formP^l We 
outline this last approach here. 

Figure [4] depicts the results of one numerical simula¬ 
tion of the initial value problem © with the solitons 
initially far apart, propagating toward each other, and 
with the shape mode unexcited, i.e I > 1, I < 0, and 
A = A = 0. Additionally A is small enough for cap¬ 
ture to occur. In this simulation, E(t) > 0 before the 
collision, so that the solution is outside the separatrix in 
Figure [2^b) , and crosses to the inside where E < 0 at 
the first collision time t\. At each subsequent collision 
E(t) jumps, reaching a plateau Ej between collisions at 
tj-i and tj, and escaping to infinity when E(t) > 0 once 
again. Upon each collision the amplitude and phase of 
A{t) jump. We represent the solution before the collision 
at time tj by the energy level Ej and by assuming 

A{t) ^ C ( cj cos u:(t — tj) Sj sinc<;(£ — tj)) (9) 
for t between tj -1 and tj , where setting the constant 

- / F(X s (t))e zut dt 
^ 4 — 00 


c r°° 

C = — F(Xs(t)) cosujtdt 

^ J — GO 




FIG. 4. Components used to construct the iterated map. 
Top X{t) showing the collision times tj. Middle: E(X, X) 
the energy of the A-component. Bottom: A{t). 


and that 


Ej+i — Ej A C 2 cd 2 - T Sj^J (11) 

The change of energy is computed using a Melnikov 
integral P*. Condition ( [lQ| ) is written in terms of (t — tj ), 
whereas the form of the map requires (t — tj+ 1 ). We 
can compute the time between collisions by matching 
between two the outer expansions connected by an in- 

_ I /Q 

ner expansion, which gives tj+% — tj = aE- + [ +o( 1) for 
some a determined by matched asymptotics. This gives, 
approximately, 


c j+1 


cos 6 j+ 1 

sin 9j+\ 


c i 



— sintfj+i 

cos6j+i 


Sj — 1 


where Qj+i = u(tj+i —tj). 

The discrete map possesses a conserved quantity % = 
Ej + j 2 C 2 (Cj A Sj) that allows us to eliminate Ej from 
the system. Defining a complex variable Zj = Cj A isj , 
the new map may be written 


allows us to scale the variables to be 0(1) in the final 
form of the system. We seek a map {Ej+ 1 , Cj+ 1 , Sj+ 1 ) = 

M{Ej , Cj ,Sj). 

The map is constructed by building a matched asymp¬ 
totic expansion which alternates between “outer expan¬ 
sions” where X(t) is approximated by As (t — tj ), Fig¬ 
ure ©b), and “inner expansions”. On the outer solution, 
valid when A (t) 1, the modes exchange energy. We 
find, by variation of parameters, that outer solutions with 
“before-collision” condition © as t — tj -A — oc, satisfy 
the “after-collision” condition 

A(t) ~ C ( Cj cos uj{t — tj) + (sj — 1) smuj(t — tj)) 

as t — tj — y Too, (10) 


^• +1 =e- i “(«-^- 1 1 2 )" (zj- 1). 


( 12 ) 


Figure [l] shows that the map reproduces the qualita¬ 
tive and many quantitative features of the ODE system 
and enables the determination of many of the features 
of the dynamics of system ©• For example, in the ex¬ 
periments depicted in Figure [l] the initial condition is 
Zj = 0. Therefore, we may estimate the critical velocity: 
the solution is captured if E\ < 0, i.e. if v < v c = ooC, 
making use of equation (11). The analysis also allows 


us to calculate approximate formulas for the centers of 
the two- and three bounce resonance windows, allowing 
the computation of the constants in equation 0 without 
data fitting.^ 





















III. A PHYSICAL MODEL OF THE COLLECTIVE 
COORDINATE EQUATIONS 
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We wish to design a surface such that a ball confined to 
roll along that surface satisfies equations of motion simi¬ 
lar to system ^ . We derive the evolution equations cor¬ 
responding for motion along a general surface and then 
design a surface that gives the desired dynamics. 


A. The mathematical model 

Consider the behavior of a point particle of mass m 
confined to a surface z = h(x,y) and moving under the 
influence of constant gravity in the ^-direction, and ignor¬ 
ing any friction or other dissipative mechanisms. We do 
not model in detail the rotation of the ball, instead not¬ 
ing that the rotational kinetic energy of a sphere rolling 
with speed v is | mv 2 . Adding this to the translational 
component gives a total kinetic energy 

T= l^> (x 2 + V 2 + z 2 ) = 5 y (^ 2 + y 2 + + hyV) 2 ^) 

The gravitational potential energy is just U = mgz = 
mgh(x,y). The evolution of a Lagrangian system with 
coordinates qj is governed by the Euler-Lagrange equa¬ 
tions 

Tt7T~7r =0 where L = T ~ U - 

at OQj OQj 

For this system, this gives, after some algebra, 


(13) 


X 

h xx x T 2h xy xy T hyyij T g 

fox 


"O' 

y_ 

\ + h 2 x + h 2 y 

foy_ 


0 


where g = Under the assumption that the height 
h(x,y ) varies slowly (e.g. h(x,y) = h(5x,Sy), { « 1), 
this is approximately 


x + gh x (x,y) = 0; 
y + gh y (x,y ) = 0. 


(14) 


Thus, when 


cJ 


g h(x, y) = U (x) + y y 2 + eF(x)y, 


(15) 


equation ( |14| ) is identical to equation Such a sur¬ 
face is shown in Figure ([5| . Numerical simulations of 
both equations (13) and ( [14] ) show qualitatively the same 
chaotic scattering behavior(not shown here). Figure[5|a) 
shows a contour plot of the surface used. 

We denote by x and y , respectively, the longitudinal 
and transverse, directions. To interpret the ODE, con¬ 
sider the case e = 0, in which case the subspace y = y = 0 
is invariant. A ball that begins rolling down the center 
of the surface with y = 0 stays along that line. When 
e > 0, this symmetry is broken, deflecting the longitudi¬ 
nal trajectory in the transverse direction. We will say a 
“bounce” occurs when the ball reaches a minimum in the 
longitudinal direction. 




FIG. 5. (a) Contour plot of the energy surface, (b) Photo¬ 
graph of the experimental surface. 


B. Effects of dissipation 


To account for experimental energy loss, we include a 
viscous damping force Fd amp = —pb(x,y), which modifies 
the equations of motion to: 

5/i . x 

x-h ——x + gh x {x, y) = 0; 

7m 

5/i . „ . . 

y + —y + gh v (x,y) = 0. 

In Goodman, Holmes, and WeinsteirP^, a dissipative cor- 
rection to system <© is derived to account for the ef¬ 
fect of radiative damping in nonlinear wave collisions. It 
adds a term like F(X) 2 A 2 A to the left-hand side of equa¬ 
tion (6b). This damping term is strongly localized in X, 


only applying when the kink and antikink are close to¬ 
gether, and nonlinear in A so that small radiation damps 
very slowly. Both types of dissipation alter the fractal 
structure displayed by the ODE system, as seen in fig¬ 
ure [6j but their effect is rather different. The localized 
damping preserves much more of the fine structure of 
the chaotic scattering, whereas the viscous damping de¬ 
stroys almost everything but the two-bounce resonance 
windows. Their behavior for v near v c is also quite dif¬ 
ferent, with viscous damping reducing the final speed of 
solutions in the two-bounce solutions much more than 
the localized damping. We therefore expect to see mostly 
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FIG. 6. Input/output plot with (a) linear damping, y — 0.01 
and (b) localized nonlinear damping, y = 0.08. 


the two-bounce resonance phenomenon in our physical ex¬ 
periments,, rather than the full chaotic scattering picture. 


Numerical experiments were used to both verify that 
system (13) could reproduce the dynamics demonstrated 
by solitary wave collisions. These produced the expected 
results, so we proceeded with laboratory experiments. 


IV. THE LABORATORY EXPERIMENT 


To test experimentally whether the physical system 
modeled by (13) displays chaotic scattering, we needed to 


fabricate a surface satisfying (15). The surface was milled 


out of high-density urethane foam using a three-axis mill 
in the Fabrication Laboratory of the NJIT Department 
of Architecture. Dimensions of the surface are given in 
Appendix [A] The rough milled surface was sanded and 
painted, as shown in Figure [5|b). A ramp was placed at 
one end, and a rubber coated steel ball (from a computer 
mouse), was rolled down this ramp and allowed to move 
along the surface until either (1) it became clear that it 
would not escape, or (2) the ball returned close enough 
to the starting point that we deemed it to have escaped. 

We computed an effective friction constant by the fol¬ 
lowing procedure. The ball was placed at the edge of 
the channel near the left edge of figure 11 and allowed 
to move under the influence of gravity. Coupling to the 
x-direction is weak here, so the motion remains largely 
confined to the ^/-direction over the time scale of observa¬ 
tion. From the video we find the sequence of successive 
maxima \y{tjf) | from which we can find the frequency and 
decay rate. Nondimensionalizing using these scales, gives 
y motion satisfying the non dimensional equation 

y + 0.019 y + y = 0. 

This nondimensionalization has a time unit of 0.15 sec¬ 


onds. The video trials shown below last on the order 
of 5-10 seconds, which is 30-60 time units, over which 
dissipation may be significant. 

In the initial experiments, we used a large marble, 
which lost energy rapidly. The rubber-coated steel ball, 
being heavier and quieter, kept rolling for much longer. 
Dissipation was not negligible, however, so any trial in 
which the ball did not exit after three “bounces” was 
considered to be trapped. 

The motion was recorded at 125 frames/second with 
a Photron FASTCAM 1024 PCI high-speed camera, po¬ 
sitioned above the surface and pointing downward. Its 
high frame rate was not necessary, but its short expo¬ 
sures prevented the ball’s image from being smeared into 
a snake-like shape. The videos were analyzed using the 
MATLAB Image Processing Toolbox, giving a time se¬ 
ries for its x and y coordinates. A screenshot of the video 
and the running program are shown in Figure [7] 

The initial velocity was varied by varying the height 
of the ball’s release. Due to the size of our foam block, 
the surface could only be milled to a fairly shallow depth. 
This meant that a ball with sufficiently large initial veloc¬ 
ity would simply fly off the milled portion of the surface 
and, unfortunately, left us unable to determine the crit¬ 
ical velocity v c experimentally. Figure [8] shows several 
interesting trajectories: first, five two-bounce trajecto¬ 
ries with the number of transverse oscillations increasing 
from two in subfigure (a) to six in subfigure (e). Subfig¬ 
ure (f) shows a fit of the time between the two bounces 
vs. the number of transverse oscillations, showing the 
same approximate linear fit as given in equation 0 , in 
this case T n « 0.961 n + 0.755. Unfortunately, the ini¬ 
tial velocities that led to these solutions do not increase 
monotonically as predicted by the theory, due, perhaps, 
to our inability to control with sufficient precision the 
transverse component of the initial location and velocity 
vectors. The figure also shows (g) a trapped trajectory 
and (h) a three-bounce trajectory. These represent, as 
far as we knox, the first (non-numerical) experimental 
evidence of the two-bounce resonance phenomenon, and 
a trace of chaotic scattering. 


V. CONCLUSIONS AND PLEA 

We have designed and built an experimental system 
that is approximately governed by the same reduced sys¬ 
tem that gives rise to the two-bounce resonance phe¬ 
nomenon in solitary wave collisions. Experiments run 
with this apparatus qualitatively reproduce the dynam¬ 
ics seen with solitary waves. 

Nonetheless, we remain hopeful that someone will find 
a way to demonstrate this phenomenon in a laboratory 
setting with actual solitary waves. The basic requirement 
is a solitary wave that supports an additional degree of 
freedom to which it can transfer energy, and which can 
transfer energy back. Usually, this takes the form of a 
mode that is localized near a stationary potential 19 22 











































FIG. 7. (Left) One frame from a movie of the experiment. (Center) The estimated trajectory up to time t « 2. (Right) The 
estimated coordinates as a function of time. See animation here. 
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FIG. 8. In all but (e) Longitudinal position x(t) (solid, left y- 
axis) and transverse coordinate y(t) (dashed, right y- axis) for 
selected initial velocities, (a-e) Two bounce trajectories. Each 
contains one transverse oscillation more than the one that 
precedes it. (f) Time between maxima of x(t) versus number 
of transverse oscillations, with best-fit line, (g) “Captured” 
trajectory, (h) Three bounce trajectory. 


or a n inte rnal mode of oscillation that moves with the 
wavep^DS although this is not always necessaryP^ 30 
Combing through some of the literature on this phe¬ 
nomenon, we find that many papers discuss the phys¬ 
ical systems described by their equations and mention 
that—perhaps—this phenomenon may be found in these 
systems. Experimental verification for some of these sys¬ 
tems is clearly impossible: Anninos et alP^l describe the 
c p A equation as describing the interactions of large-scale 
domain walls in the universe. They admit, “Because the 
possibility of head-on collisions is small in the real Uni¬ 
verse, our findings are not expected to have a profound 
effect cosmologically.” At the other extreme of scales, 
Kudryavtsev describes the ip 4 system as a model for the 
Higgs field. 1031 Other unpromising applications for the 
theory are information transport in brain microtubules 
and various other quantum field theories! 32 * 33 * 

Other applications are perhaps closer to being exper¬ 
imentally realizable. In their original paper on the sub¬ 
ject, Campbell et al. suggest excitations in polymeric 
chains and phase transitions in uniaxial ferroelectrics! 13 ^ 
Tan and Yang study the phenomenon in vector solitons 
in optical fiber e t al.’s numerical experiments 
potentially describe the scattering of solitons off a de¬ 
fect in a long Josephson junction. 19 22 Finally, Forinash 
et al. mention the denaturation of DNA.^This list is far 
from exhaustive, and we encourage readers to think if an 
experiment is possible in their favorite system. 
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Appendix A: Specification of the experimental apparatus 

To design the experimental surface, we performed nu¬ 
merical experiments of system (Ml), varying not just the 
parameters in the potential d 6 ), but also the formulas 
for the potentials U{x) and F(x), trying to ensure we 
could observe interesting dynamics using the materials 
available to us. The surface chosen was 

h(x, y) = rj(e ax - e bx ) + cy 1 2 3 + eye dx 

where the units of distance in all coordinate directions is 
centimeters. The parameters chosen were 

a = 0.389 cm -1 , b = 0.306 cm -1 , c = 0.0764 cm -1 , 
d = 0.306 cm -1 , 77 = 3.27 cm, e = 0.25. 
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